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Abstract 

We consider an idealized model in which individuals' changing opinions and their social network 
coevolve, with disagreements between neighbors in the network resolved either through one imitat- 
ing the opinion of the other or by reassignment of the discordant edge. Specifically, an interaction 
between x and one of its neighbors y leads to x imitating y with probability (1 — a) and otherwise 
(i.e., with probability a) x cutting its tie to y in order to instead connect to a randomly chosen 
individual. Building on previous work about the two-opinion case, we study the multiple-opinion 
situation, finding that the model has infinitely many phase transitions. Moreover, the formulas 
describing the end states of these processes are remarkably simple when expressed as a function of 
(3 = a/(l — a). 
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I. INTRODUCTION 



In the last decade, there have been a number of studies of systems in which the states of 



individuals and the connections between them coevolve, see [l 



the spread of opinions 




include evolutionary games [3||-[7f and epidemics |s]-[l2j], but here we wil 



2j. The systems considered 



13(-[l7|. Different from the models of cascades 18j- 20[ which are 



concentrate on 



also widely used in the study of opinion spread, the evolving voter model we study here 
allows an agent to switch between different opinions and the network topology to change 
accordingly, yet we assume that agents impose equal influence over each other (cf., multi- 
state complex contagions 2JJ-23J). This model provides building blocks to quantitatively 
study collective behaviors in various social systems, e.g., segregation of a population into 
two or more communities with different political opinions, religious beliefs, cultural traits, 
etc. 



We are particular 



y interested here in systems that generalize the model proposed by 



Holme and Newman 24|. In their model there is a network of N vertices and M edges. The 



individual at vertex v has an opinion £(v) from a set of G possible opinions and the number 
of people per opinion 77V = N/ G stays bounded as N gets large. On each step of the process, 
a vertex x is picked at random. If its degree d(x) equals 0, nothing happens. If d(x) > 0, (i) 
then with probability 1 — a a random neighbor y of x is selected and we set = £(y); (ii) 
otherwise (i.e., with probability a) an edge attached to vertex x is selected and the other 
end of that edge is moved to a vertex chosen at random from those with opinion £(x). This 
process continues until the 'consensus time' r, at which there are no longer any discordant 
edges — that is, there are no edges connecting individuals with different opinions. 

For a = 1, only rewiring steps occur, so once all of the M edges have been touched, the 
graph has been disconnected into G components, each consisting of individuals who share 
the same opinion. Since none of the opinions have changed, the components of the final 
graph are all small (i.e., their sizes are Poisson with mean 7^). By classical results for the 
coupon collector's problem, this requires ~ MlogM updates, see e.g., page 57 in 25| . In 
the case of sparse graphs we consider here M ~ cN (i.e., M/N — > c) so the number of steps 
is O(NlogN), i.e., when N is large it will be ~ CN log N. 

In contrast, for a = this system reduces to the voter model on a static graph. If we 
suppose that the initial graph is an Erdos-Renyi random graph in which each vertex has 



average degree A > 1, then (see e.g., Chapter 2 of [26j) there is a "giant component" that 
contains a positive fraction, /i, of the vertices and the second largest component is small 
having only 0(log N) vertices. The voter model on the giant component will reach consensus 
in 0(N 2 ) steps (see, e.g., Section 6.9 of 26J), so the end result is that one opinion has fiN 
followers while all of the other groups are small. 

Using simulation and finite size scaling, Holme and Newman showed that there is a critical 
value a c so that for a > a c all of the opinions have a small number of followers at the end of 
the process, while for a < a c "a giant community of like-minded individuals forms." When 
the average degree A = 2M/N = 4 and the number of individuals per opinion 7^ — > 10, this 
transition occurs at a c ~ 0.46. See [27]- 30] for recent work on this model. 



In 



3l| . we studied a two-opinion version of this model in which on each step an edge 



is chosen at random and is given a random orientation, (x,y). If the individuals at the 
two ends have the same opinion nothing happens. If they differ, then (i) with probability 
1 — a we set £(x) = £(y); (ii) otherwise (i.e., with probability a) x breaks its edge to y 
and reconnects to (a) a vertex chosen at random from those with opinion £(x), a process 
we label 'rewire-to-same', or (b) at random from the graph, a process we label 'rewire-to- 
random'. Here, we will concentrate on the second rewiring option, rewire-to-random. While 
this process may be less intuitive than the rewire-to-same version, it has a more interesting 
phrase-transition, as documented in 3l| . 

The remainder of this paper is organized as follows. In Section [Til we recall the main 
results from 3l| that provide essential context for our observations of the multiple-opinion 
case, which we begin to explore in Section IIIII We then continue in Section [IV] with fur- 
ther quantitative details about the phase transitions and their underlying quasi-stationary 
distributions, before concluding comments in Section IVl 



II. TWO-OPINION MODEL 



Suppose, for concreteness, that the initial social network is an Erdos-Renyi random graph 
in which each individual has average degree A > 1, and that vertices are assigned opinions 
1 and independently with probabilities u and 1 — u. Simulations suggest that the system 
has the following 

Phase transition. For each initial density u < 1/2 there is a critical value a c (u) so that 
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FIG. 1. Simulation results for the rewire-to-random model, starting from Erdos-Renyi graphs with 



N=100,000 nodes and average degree A = 4 



311 ] . The final fractions of the minority below phase 



transitions follow an universal curve independent of the initial fractions. 

for a > a c (u), consensus occurs after 0(N log N) updates and the fraction of voters with 
the minority opinion at the end is p(a,u) « u. For a < a c (u) consensus is slow, requiring 
0(N 2 ) updates, and p(a,u) ~ p(a,0.5). 

To help understand the last statement, the reader should consult the picture in Figured! 
If the initial fraction of l's u = 1/2 then as a decreases from 1, the ending density p(a, 1/2) 
stays constant at 1/2 until a = a c (l/2) and then decreases to a value close to at a = 0. 
For convenience, we call the graph of p(a, 1/2) for a < a c = 0.74, the universal curve. If the 
initial density is u < 1/2, then the ending density p(a,u) stays constant at u until the flat 
line (o?, u) hits the universal curve and then p(a, u) ~ p(a, 0.5) for a < a c (u). The main aim 



of 31] was to use simulations, heuristic arguments, and approximate models to explain the 
presence and properties of this universal curve describing the consensus states that result 
from the slow- consensus process. To make it easier to compare the results here with the 
previous paper, we rescale time so that times between updating steps are exponential with 
rate M, where M is the total number of edges. 

Quasi-stationary distributions. Let v(a) = p(ct, 0.5). If a < a c (l/2) and v(a) < 
u < 1/2 then starting from product measure with density u of 1 's, the evolving voter model 
converges rapidly to a quasi- stationary distribution v a ^ u . At time tM the evolving voter model 
looks locally like v a ,e(t) where the density changes according to a generalized Wright-Fisher 
diffusion process 

dB t = y/(l-a)[cJ t (l-e t )-b Q ]dB t (1) 



until 6 t reaches v(a) or 1 — v{a), the two solutions of c a x(l — x) = b a . 

To further explain the phrase "quasi- stationary distributions" in this context, we refer 
the reader to Figure [2j Let Ni(t) be the number of vertices in state 1 at time t, A*oi(£) be the 
number of 0-1 edges (that is, the number of edges connecting nodes x and y with = 0, 
£(?/) = 1). Similarly, let N a b c (t) be the number of connected triples x-y-z with £(x) = a, 
£(y) = b, and £(z) = c. The top panel of Figure |2] plots N i(t)/M versus Ni(t)/N for five 
different simulations (with different initial densities, u) for a = 0.5. Note that in each case 
the simulation rapidly approaches a curve ~ 1.710a;(l — x) — 0.188 and then diffuses along 
the curve until consensus is reached (TVoi = 0). At both of the possible consensus points on 
the curve, the fraction of the minority opinion is ~ 0.12, in accordance with the simulation 
in Figure [U 

The bottom panel of Figure [2] similarly plots N 010 (t)/N versus Ni(t)/N for a = 0.5 
and u = 1/2. Again the simulation rapidly approaches a curve (approximately cubic) and 
diffuses along it until consensus is reached. Since iVoio = if iVoi = 0, and it is very unlikely 
that all 0-1's only occur in 0-1-1 triples, the zeros of the cubic curve for 0-1-0 and quadratic 
curve for 0-1 coincide. 

One can repeat the simulations in Figure [2] for other network measurements, with the 
result that their values are similarly determined by the density u(t) = Ni(t) /N . This is some- 
what analogous to a stationary distribution from equilibrium statistical mechanics — e.g., the 
Maxwell-Boltzmann distribution associating the velocity distribution with the temperature. 
We call our distributions quasi-stationary because our system is a finite state Markov chain, 
which will eventually reach one of its many absorbing states N 01 = 0, and hence there is no 
true stationary distribution. Nevertheless, an improved understanding of the system is ob- 
tained from these observations, displaying a fast dynamics rapidly converging to a family of 
neutrally-stable quasi-stationary distributions followed by slow, diffusive dynamics through 
the space local to the quasi-stationary distributions until consensus is reached. 

To begin to explain the behavior of 9t given in ([I]), note that when an edge is picked with 
two endpoints that differ, a rewiring will not change the number of l's, while a voting event, 
which occurs with probability (1 — a), will result in an increase or decrease of the number 
of l's with equal probability. When 9 t = u the rate at which 0-1 edges are chosen is equal 
to the expected fraction of 0-1 edges under i/ a>u , which is c a u{l — u) — b a . 



As shown in 



3l| , the behaviors for the rewire-to-same model in terms of quasi-stationary 
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FIG. 2. (Top) Evolution of the fraction of edges that are discordant 0-1 edges, Ani (i)/M, versus the 
population of opinions Ni(t)/N when a = 0.5 for the rewire-to-random dynamic. Five simulations 
starting from u=0.2, 0.35, 0.5, 0.65, and 0.8 are plotted in different colors. Each simulation 
starts from an Erdds-Renyi graph with N=100,000 nodes and average degree A = 4. After initial 
transients, the fraction of discordant edges behaves as a function of the population of opinions. 
(Bottom) Similarly, the number of 0-1-0 connected triples behaves as a function of Ni/N after an 
initial transient (one simulation). 

distributions are very similar, but with small differences from the rewire-to-random model 
that yield fundamentally different consensus states. In rewire-to-same, there are quasi- 
stationary distributions v' au under which the expected fraction of 0-1 edges is c' a u(l — u). 
Again the simulation comes rapidly to this curve and diffuses along it until consensus is 
reached. That is, unlike Figure |2] (Top), the arches of quasi-stationary Nqi/M values versus 
Ni/N maintain their zeros at Ni/N = {0, 1}. Thus, for a < a' c (l/2), the minority fraction 
obtained at the consensus time is always ~ for rewire-to-same. 
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III. MULTI-OPINION MODELS 



Bohme and Gross [32| have studied the three-opinion version of the evolving voter model 
with rewire-to-same dynamics. In this case, the limiting behavior is complicated - one may 
have partial fragmentation (l's split off rapidly from the 2's and 3's) in addition to full 
fragmentation and coexistence of the three opinions. See their Figures 3-5. As we will see 
in the present section, the behavior of the multi-opinion rewire-to-random model is much 
simpler because small groups of individuals with the same opinion will be drawn back into 
the giant component. We thus aim to extend the understanding of the two-opinion model 
behavior to larger numbers of opinions. 

Consider now the /c-opinion model in which voters are assigned independent initial opin- 
ions that are equal to i with probability ttj. Let u = (u±, itfc) and let N-l be the 
number of edges at which the endpoint opinions differ. When k — 3, frequencies of the three 
types must lie in the triangle of possible values A = {u = (ui,u 2 ,u 3 ) : Ui > 0, Yli u i = !}• 
To preserve symmetry, we draw A as an equilateral triangle in barycentric coordinates by 
mapping (x,y, z) — > (x, z\/3/2). The top panel in Figure [3] plots N^(t)/M as a function of 
the opinion densities as the system evolves, generalizing the one- dimensional arch observed 
for k = 2 to a two-dimensional cap for k = 3. 

Generalizing the parabolic form of the arch for k = 2, we conjecture 

E U N^/M = I 



As in the two opinion case, the simulated values come quickly to the surface and then diffuse 
along it. In some situations, one opinion is lost before consensus occurs and the evolution 
reduces to that for the two opinion case. However, in one of the simulations shown, the 
realization ending with x ~ 0.5, all three opinions persist until the end. 

The picture is somewhat easier to understand if we look at the cap from a top view, where 
the E u N^l = level sets for different a are observed to be circles. In the bottom panel of 
Figure |3] we plot the E U N^ = circles for different a's fitted from simulation data using Eq. 
(J2J) as well as the consensus opinion frequencies from the simulations (indicated by small 
circle data points). The two agree with each other up to small stochastic fluctuations. The 
size of the E U N^ = level set then dictates different consensus state properties. For example, 
the circle corresponding to a = 0.5 intersects A in three disconnected arcs. As a increases, 
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FIG. 3. (Color Online) Top: plot of the fraction of discordant edges versus the population of 
opinions in barycentric coordinates for three opinions and a = 0.5. Multiple simulations corre- 
sponding to different initial densities are shown while each one starts from an Erdos-Renyi graph 
with N=10,000 nodes and average degree A = 4. Similar to the two-opinion case, the simulations 
quickly converge to a parabolic cap of quasi-stationary distributions. Bottom: top view of the 
parabolic caps of quasi-stationary distributions for q=0.1,0.2,...,0.8. We fit the parabolic cap Eq. 
([2]) to simulation data at various a's and then plot the level sets E U N^ = 0, which are the inter- 
sections of the parabolic caps with the = plane, as the large circles with colors indicating 
values of a. 

the radius of the E U N^ = level set decreases. When a > a c (l/2), the critical value of the 
two opinion model, the circle E U N^ = falls fully inside the triangle, so an initial condition 
including all three opinions will continue to demonstrate all three opinions at consensus. 
For example, the small circles around the innermost circle give the ending frequencies for 
several simulations for a = 0.8. If the initial frequencies fall within the E U N^ = circle, 
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then the model will quickly relax to the quasi-stationary distributions above the circle and 
then diffuse along the cap until consensus is reached at some E u N-t = point. If instead 
the initial frequencies u fall outside the E U N^ = circle — that is, for a above the phase 
transition point a^(u) — the consensus time jumps from 0(N 2 ) to O(NlogN), similar to 
a c (u) for the two-opinion model, with the final opinion frequencies essentially the same 
as the initial u. What is new in this case is that when starting with three opinions and 
a c {u) < a < ctziu) < a3({|, |, |}), the system always ends up with three distinct opinions. 

For k > 3, our simulation results indicate the same type of behavior as the system evolves. 
We define a& to be the largest a for which consensus takes 0(N 2 ) updates when we start 
with k opinions with density 1/k for each opinion. Then as k — > oo the multi-opinion 
model has infinitely many phase transitions. When < a < consensus occurs after 

0(N log N) steps if we start with k opinions, while if we start with k + 1 equally likely 
opinions the system quickly converges to a quasi-stationary distribution and diffuses until 
consensus occurs after 0(N 2 ) updates and there will always be k+ 1 opinions present at the 
end. The associated picture is the natural dimensional extension of the relationship between 
the k = 2 and k = 3 models: just as a.2 = a c (l/2) corresponds to the point at which the 
E u N-l = circle for k = 3 is the inscribed circle within the A triangle, a-s corresponds to 
the point at which the E u N-l = circle reaches zero radius — that is, the point at which 
the E U N^ = sphere for k = 4 has become the inscribed sphere within the corresponding 
barycentric tetrahedron. 



IV. QUANTITATIVE CHARACTERIZATION OF QUASI-STATIONARY DIS- 
TRIBUTIONS 

For each k we simulate our multi-opinion rewire-to-random model starting from k opin- 
ions with each opinion taking 1/k fraction of nodes at random for a wide range of a's. 
Generalizing the picture of the one- dimensional arch for k = 2 and the two-dimensional 
cap for k = 3, the number of discordant edg function of frequencies conjectured in 

Eq. ([2]) is a co-dimension 1 hypersurface characterizing the quasi-stationary states, and the 
behavior of the equal-initial-populations case will allow us to describe this surface, thereby 
characterizing behaviors for general initial populations. 

First the critical a^'s are identified when the slow diffusion of iV^ cannot be observed 
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FIG. 4. Coefficient cq{0) (left) and C2(/3) (right) in Eq. ([2]) for models with multiple opinions. 
Each value of the coefficients is obtained by fitting Eq. ([2]) to multiple simulations starting from 
Erdos-Renyi graphs with N=100,000 nodes and average degree 4. The fitting error is very small 
(R 2 « 0.99) except for close to the critical values. 

for the first time as a increases from to 1. Then we fit N^(t)/M to Ui(t) = Ni(t)/N 
(i = l,...,k) using Eq. (j2]) at every a up to a^, and plot the fitted coefficients Co and c 2 
against = a /{I — a) in Figure HI Remarkably, the coefficients in (j2J) appear to be well 
approximated by linear functions oi = a/(l — a). The graphs shows some curvature near 
— 0, which may be caused by the fact that — (a = 0) corresponds to a voter model 
without evolution of the underlying network. In the rest of the paper, we will work with 
for simplicity. Naturally, critical points translate to 0k = «fc/(l — «&)• 

The fitted coefficients from the 2-opinion model deviate slightly from those fitted from 
higher-order models, which implies that Eq. ()2]) is not universal for the multi-opinion model 
and higher-order terms are possible. However, while the discrepancy between the fitted 
coefficients of the 2-opinion model and those of the 3-opinion one is apparent, difference be- 
tween fitted coefficients of higher-order models is negligible, which implies that the inclusion 
of higher-order terms beyond the 3rd would not make significant changes to the equation. 
To probe the effect of higher-order terms we introduce terms up to kth order for k opinions. 
Noting (J2i u i) 2 = 1) Eq. © is equivalent to: 

k 

(3) 



E U N^/M = -c {a) + c 2 {a) 



U{Uj . 



Given the symmetry of the system in itj's, the only possible choice in degree-k polynomials 
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FIG. 5. Coefficients cq{(3) (left) and c 2 (/3) (right) in Eq. for models with multiple opinions. 
Each value of the coefficients is obtained by fitting Eq. @ to the same data as in Figure HI 



is: 



E U N^/M = -c (a) + c 2 (a) ^ u h i 

{h,i2}eA2 

+c 3 (a) 22 u ii u h u i3 

+c fc (a) ^ u h u h - ■ -u ik , 



(4) 



where *4.j is the collection of all z-element subsets of {1, 2, /c}. Using the same simulation 
data as above, we refit N^(t)/M to Wj(i)'s (z = 1, k) according to the generalized formula 
Eq. (J2J) and plot the fitted coefficients c and c 2 against /3 in Figure Fitting diagnostics 
suggest that higher-order terms are significant from zero (with p- value < 10 -4 ) and it can 
be seen that those terms explain the inconsistency between fitted coefficients of different 
models in Figures HI However, the difference between the two fitted functions of Eq. ([2]) and 
Eq. (j4]) is actually small (w .1 in L 2 -norm) and thus higher-order terms are small corrections 
to the hyper-surface Eq. (j2J). 

Values of the coefficients q(/3) for the three opinion model near its critical value f3% ~ 5.2 
show some scatter, but this is to be expected since the surface is very small at this point. 
Values for the four opinion model appear to become more difficult to fit prior to /?4 since 
E u N-l = is a three-dimensional hyper-surface in four- dimensional space, so much more 
data is required to get reliable estimates of coefficients. 

As is visually apparent in Figure the coefficients c and c 2 for the first two terms 
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FIG. 6. Coefficient cs(f3) in Eq. ([4]) for models with multiple opinions. Each value of c^{f3) is 
obtained by fitting Eq. @ to the same data in Figure HI 

in Eq. (j4j) are well approximated by linear functions, with best fits Co(/3) ~ 0.22/3 and 
c 2(/3) ~ 1.3 + 0.38/3, while coefficients for higher-order terms are not linear in /3 (e.g., see 
Figure [6] for c${f3)). For comparison, the best fits for cq and c 2 in Eq. (j2J) (as in Figure H]) 

are 

c 2 (/3) « 1.3 + 0.5/3, c (/3) « 0.25/3. (5) 

Since Eq. (j2J) well approximate the higher-order hyper-surface Eq. (j3j), its simple form can 
be used to estimate the critical points for phase transitions. Combining (j2j) and (111), and 
then solving 

(0.65 + 0.25/3)(l - k(l/k) 2 ) - 0.25/3 = 

gives 

fa = 2.6(fc - 1) . 

which agrees with the critical ^'s identified when the slow diffusion of cannot be observed 
in simulations as increases. 

V. CONCLUSION 

Our multi-opinion voter model has infinitely many phase transitions. When /3& < /3 < 
Pk+i, consensus occurs rapidly when we start with k opinions, while if we start with k + 1 
equally likely opinions there will always be k + 1 opinions present at the end. To a good 
approximation /3 fc = 2.6(k — 1), but the departures from linearity in the plots of c 2 (/3) 
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and co(/3) suggest that this result is not exact. However, formulas for various quantities 
associated with this model are close to polynomials, so an exact solution may be possible. 

More complicated rewiring rules might also be considered, particularly if they main- 
tained high clustering or other global macroscopic properties. An even more complete un- 
derstanding of the present rewiring system would help motivate similar investigations for 
other rewiring rules. 
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